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ON THE GROUND-STATE WAVE FUNCTION OF THE ONE-DIMENSIONAL 
POLARON IN THE STRONG-COUPLING LIMIT 

ROHAN GHANTA 


Abstract. We consider the one-dimensional Frohlich polaron localized in a symmetric decreasing electric 
potential. It is known that the non-linear Pekar functional corresponding to our model admits a unique 
minimizer. In the strong-coupling limit, we show that any approximate ground-state wave function of our 
model- after integrating out its phonon modes- converges in the weak sense to this unique minimizer. 


1. Introduction 

As a model of an electron moving in an ionic crystal, the polaron continues to be of interest. Because it 
is also one of the simplest examples of a particle interacting with a quantum field, it has served as a testing 
ground for various techniques in field theory (see [AlDe]) such as the Feynman path integral (see [Fy]). But 
despite the attention it has received over the last eight decades, many questions remain open. With a few 
exceptions in a limiting case (see e.g. [DoVa], [LiTh]), the polaron has eluded the exact calculation of the 
most basic quantitites such as the effective mass and the ground state energy. Moreover, an exact analytic 
expression for the ground-state wave function of the model has yet to be given. In this paper, we give the 
first convergence result for the wave function. 

Polaron theory began in the 1930s when ionic crystals of nonmetallic type- capable of producing strong 
electric fields- were being introduced into electronic devices (see [Pek]). Starting with L.D. Landau’s paper 
[Ld] from 1933, where it was suggested that the electron deforms the crystal and traps itself in a hole of its 
own making, various models (see [Dev]) were developed to explain the transport of electrons through these 
crystals. These models were considerably more reliable for experiments than the standard band theory (see 
[Pek]), because they account for the polarization of the crystal due to the moving electron. This polarization 
is modelled in terms of the vibrational displacement of the ions in the crystal lattice; in the literature these 
vibrations are called phonon modes. 

In 1937 H. Frohlich suggested a model- known today as the Frohlich polaron- to explain electrical break¬ 
down in these crystals [Fr]. The Hamiltonian is 


k 



E 


^ ikx 

|k| 



( 1 . 1 ) 


where p = —iV (the electron momentum) and acts on -F(8)L^(]R^), where J- denotes the (symmetric) phonon 
Fock space. In (11.11) we use “cc” to denote the electronic coordinate, “fc” for the phonon mode and F for the 
volume of the crystal. The creation and annihlation operators a\ and Ofc are defined on J- with the canonical 
commutator relation [a^, a\,\ = 6{k — k'). The coupling paramter a > 0 was introduced by Frohlich [Fr-2] in 
1954 to describe the interaction between the electron and the phonon modes. The model is also known as the 
large polaron, because the spatial extension of the wave function is larger than the crystal lattice spacing. 
Therefore a continuum approximation > r(27r)“^ f d^k for the Hamiltonian in (ED is also allowed. 

A proof of the self-adjointness for Hamiltonians of this type was first given in 1964 by E. Nelson [Ne]. 

The ground state energy of the model is 


Ea — inf {(di, iJctd')jr|gii2(R3) [ d* e .?’(8>L^(M^) and ||d'||jFg,i2(R3) — 1} , (1.2) 


where the operator Ha is the Hamiltonian in (11.11) . A normalized wave function in IF that achieves 
the ground state energy in ()1.2I) is a ground-state wave function. 
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The mathematical difficulty in calculating the ground state energy and the ground-state wave function 
of the model stems from the electron-phonon interaction term of the Hamiltonian in (HD. Not only are 
the electron and phonon co-ordinates in HD inseperable, but we also do not a priori know the explicit 
dependence between the co-ordinates. The minimization problem in (II.2p is therefore intractable. 

During the 1950s this mathematical difficulty motivated physicists to develop various techniques for 
approximating the ground state energy in (HD by exploiting the properties of the ground-state wave function. 
In his 1951 monograph [Pek] S.l. Pekar presents a non-linear theory posited on his Produkt-Ansatz for the 
ground-state wave function. Based on his physical intuition that the phonons are not sensitive to the 
instantaneous position of the electron, Pekar proposed that the ground-state wave function in (11.211 can be 
expressed in the product form 


'P = I</')®IC)- (1-3) 

In (11.31) 1C) £ -F is a coherent state defined only on the phonon co-ordinates, and (j) £ L^(K^) is a normalized, 
electronic wave function. In particular, (^) is a minimizer of the non-linear problem: 

ep = inf|/’ |VC)pdx- /" f I f ^(x)2dx = l'l. (1.4) 

J iR3xR3 1^ “ y| Jr 3 j 

Pekar’s ansatz in (II. 3p offers the computational advantage of eliminating all of the phonon co-ordinates in the 
optimization problem from (11.21) for the ground state energy. Minimizing (ik, HaA^) over the more restrictive 
set of product wave functions in (11.31) yields the following upper bound for the true ground state energy: 

Ea < inf I ||4'||jrg,p2(R3) = 1 and 4' = \4)) ® |C)} (1.5) 

= a^ep, ( 1 . 6 ) 

with the Pekar energy ep as defined in (11.41) . Numerical work in 1976 by S.J. Miyake suggests that ep = 
-0.108513 [My]. 

Pekar’s ansatz gives rise to a widely studied minimization problem with a non-linear energy functional, 
given in (HD above, known as the Pekar functional (see eg. [GrHtWl], [Li], [Ln[, [Ln-2[ and [LwRg[). That a 
minimizer actually exists for the non-linear problem in (11.41) has been shown only in 1977 by E.H. Lieb using 
rearrangment inequalities [Li[. Lieb has also established that this minimizer is unique up to translations by 
proving that there is a unique solution for the corresponding Euler-Lagrange equation: 

{-A - 2^^ \(j){y)\'^\x - y\-^dy^ (t){x) = Mix). (1.7) 

The equation in (HD is known as the Choquard-Pekar or the Schrodinger-Newton equation, and it has 
attracted a lot of attention in the recent literature (see eg. [GrnVs], [ LwRg[, [Lz[, [MzVs[, [MzVs-2[, [MzVs- 
3[, [MzVs-4[, [MzVs-5[ and [Rd]). In particular, E. Lenzmann has shown that around the unique minimizer 
Q(x) of the Ghoquard-Pekar equation in (11.71) . the linearization 

L+C = -AC + AC - {\x\-^ * m C - 2Q (|x|-i * (QC)) 

has a non-degenerate kernel [Lz[: 

kerL+ = spaiijS^iQ, da;^Q}. (1.8) 

With this non-degeneracy result, Lenzmann establishes the uniqueness up to translations of the pseudo- 
relativistic version of the Ghoquard-Pekar equation via an implicit function-type argument [Lz[. Lenzmann’s 
non-degeneracy result in (II. 8|) has found many applications (see [Rd] and the references therein), and it is a 
crucial ingredient, for example, in the proof of the symmetry of the bipolaron bound state given by Frank, 
Lieb and Seiringer in [FrLiSr]. These uniqueness and non-degeneracy results have been extended to the 
anisotropic polaron (see [LwRg]) by J. Ricaud in recent work [Rd]. 

Despite the mathematical interest that the Pekar functional continues to stimulate, it was noticed already 
in the 1950s that Pekar’s ansatz is only physically sensible for describing crystals with a very large (“strong”-) 
coupling paramter a [Fr-2]. And a more adequate theory was needed to explain the transport of electrons 
in weak-coupling crystals (eg. InSb; cf. [Dev]). In 1954 Frohlich has introduced a weak-coupling theory (see 
[Fr-2]) based on the canonical transformation of T.D. Lee, F.E. Low and D. Pines from 1953 [LLP]. This in 
turn motivated R.P. Feynman in 1955 to use the path integral to develop an intermediate-coupling theory 
that is applicable to a wider range of coupling paramters [Fy]. 
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These theories- each useful at different values of the coupling parameter- only provide an upper bound for 
the true ground state energy in (11.21) . In 1958, however, E.H. Lieb and K. Yamazaki arrive at a rigorous 
lower bound for the ground state energy via modifying the Hamiltonian dm) rather than making an ansatz 
about the ground-state wave function [LiYa]. Their lower bound, however, differs from Pekar’s upper bound 
in (11.61) by a factor of 3, so their calculation unfortunately does not yield the exact ground state energy in the 
strong-coupling limit. Lieb and Yamazaki’s techniques from 1958 nevertheless continue to inspire rigorous 
calculations in strong-coupling theory: commutator estimates with their vector operator Z = (^ 1 ,^ 2 , .^ 3 ), 
where 



ikx 

\k\>K ' ' 


J = 1,2,3 


are necessary when using an ultraviolet cutoff on Hamiltonians such as that in (11.11) (see eg. [AnLa], [BeBl], 
[FrGs], [FrSl] and [LiTh]). It is also useful in quantum electrodynamics (see proof of Corollary 2.2 in [LiLo-2]). 

To this day it is not known how to calculate the exact ground state energy at finite values of the coupling 
parameter a, but in 1981 M.D. Donsker and S.R.S. Varadhan have succeeded in using techniques from large 
deviation theory to show that Pekar’s approximation 1 ()1.5I) . ()1.6I) 1 of the ground state energy is exact in the 
strong-coupling limit [DoVa]: 

lim ^ = ep. (1.9) 

Of course, ep is the upper bound from (11.41) that is calculated using Pekar’s ansatz for the ground-state 
wave function. In 1995 F.H. Lieb and L.F. Thomas provide an alternate, pedestrian proof of dm using 
simple modifications of the Hamiltonian, a philosophy that can be traced back to the earlier work [LiYa] 
of Lieb and Yamazaki from 1958. In [LiTh] Lieb and Thomas use coherent states to obtain an agreeable 
lower bound, and the strategy also yields a rate of convergence for the result in (11.91) . Their technique has 
motivated recent study of the ground state energy of multi-polaron systems (see [AnLa],[BeBl] and [GrMl]) 
and can also be adapted to other models (see [FrGs] for the polaron in a large magnetic field). We use their 
methods to argue the exact ground state energy of our one-dimensional model (see Theorem 1 below). 


1.1. Motivation. Showing Pekar’s product wave function yields the exact ground state energy in the strong¬ 
coupling limit is a successful chapter in polaron theory, spanning more than four decades. This by no means 
is a justification of Pekar’s ansatz for the ground-state wave function, which, as Lieb and Yamazaki point out 
in 1958, is in fact inadequate for calculating the expectation values of various operators at the ground state. 
For example, the expectation value {H^) = 00 when one takes the ground-state wave function to be Pekar’s 
product function in ()1.3I) . It still remains to understand the connection between the true ground-state wave 
function and Pekar’s product wave function. 

With this paper we present a strategy that can be used to show that in the limit a —00 the true ground- 
state wave function of the Frohlich Hamiltonian- after integrating out its phonon modes- converges (in the 
weak sense) to the electronic wave function in Pekar’s ansatz. As evident from m, this electronic wave 
function is the minimizer of the corresponding non-linear Pekar functional. 

A much stronger relationship between the ground state of the Frohlich Hamiltonian and the product wave 
function from Pekar’s ansatz in (11.31) has been recently conjectured by E.H. Lieb and R. Seiringer [LiSr]. 

To have a well-defined notion of convergence, however, we need to ensure that the corresponding Pekar 
functional admits a unique minimizer. Note that the Pekar energy functional in (11.41) . originally studied by 
Lieb in [Li], has translational invariance. To even have a chance at uniqueness we must break translation 
invariance by introducing a localizing electric potential (in the x co-ordinate only) for the system: for some 
V (x) > 0 that vanishes at infinity, consider the localized Hamiltonian 

Ha{V)= He-a^V{ax), ( 1 . 10 ) 

where is the translation invariant Hamiltonian from HID. The corresponding Pekar energy functional is 

£v{(l>)=[ dx- f [ l‘?^(y)l dxdy. (1.11) 

Js.3 \ y j yR3xR3 — y| 

In 1977 E.H. Lieb proved uniqueness up to translations [Li] for the functional in (11.41) (see also “Appendix 
A” in [Lz] for a somewhat different proof), but his proof exploits some essential scaling relations that are 
no longer true for the functional in ([l.llL where we break translational invariance of the system with a 
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localizing potential. The existence of a minimizer for the Pekar functional with the inclusion of an external 
potential, (11.111) . has been studied by P.L. Lions in the 1980s when developing his compactness arguments 
(see [Ln] and [Ln-2] and cf. [Li-2], an earlier similar result of E.H. Lieb). In the presence of a localizing 
potential, while existence issues can now be settled using Lions’ concentration compactness lemma [Ln], we 
are not aware of any tools for addressing the uniqueness of the minimizer. Even when the localizing potential 
in (ll.ll|) is the Coulomb potential, |a;|“^, we are unable to adapt the technique from [Li]. 

A uniqueness result is available only in one-dimension (see [JjSt] and [MlStTr]), so we can only provide 
rigorous arguments for a one-dimensional model. But we emphasize that there is nothing intrinsically one¬ 
dimensional about our strategy, which follows from a very simple application of the variational principle. 


2. The One-Dimensional Model and Statement of Results 

We work with the one-dimensional Frohlich polaron localized in a symmetric decreasing, C'^(M) electric 
potential V (x) > 0 that vanishes at infinity. The Hamiltonian for our model. 


i?a(R) = 




+ Y. - 

| fe |>0 


i: [ 




Akx 


+ ale 


t ^—ikx 


| fe |>0 


a^V{ax), 


( 2 . 1 ) 


is defined on O L^(M), where is a symmetric Fock space over {hjL). In (12.11) we use “x” to denote 
the electronic coordinate, “k” for the phonon mode, “2L” is the length of the crystal and “i” is the lattice 
spacing. The creation and annihlation operators a\. and ak are defined on T with the canonical commutator 
relation [ak, a\i\ = 6{k — k'). As usual the continuum approximation ^ f allowed. 

The ground state energy of the model is defined as 


E^{V) = inf{(vl-, iJ„(P)vIr)^^i 2 I = 1}, (2.2) 

and an optimizing function in (12.21) is the ground state wave function of the model. As with the three- 
dimensional case discussed above, it is straightforward to calculate with Pekar’s Produkt-Ansatz (see m 
above) for the ground-state wave function 

'I' = |u)l2(r) (g) |^)jF, (2.3) 

that 

E^V) < a^e(V). (2.4) 

A scaled version of the electronic wave function in Pekar’s ansatz, is a minimizer of the one¬ 

dimensional minimization problem for the Pekar energy e(P): 


e(P) = inf{fy(u) : 

with the one-dimensional Pekar functional: 

/ v?dx = 1}, 

/R 

(2.5) 

£v{u) = f {u^ — 

— P(x)u^) dx. 

(2.6) 


R 


When P = 0, it has been shown that e(0) = —-E (see e.g. [Gh]). Since P is a localizing potential, e(P) < 0. 

The one-dimensional polaron was introduced by E.P. Gross in his 1976 paper as a toy model [Go]. But it 
has since attracted sizeable attention in the literature and is not entirely artificial (see eg. [FrGs], [FrLiSrTh], 
[Gn], [KoLeSm], [Mn], [PtSm], [SKVPD], [Sp], [Sp-2], [VPSD]). The most compelling reason to study the 
one-dimensional model is its connection to the three-dimensional polaron in a magnetic field, also known 
as the magneto-polaron. In 1992 E.A. Kochetov, H. Leschke and M.A. Smondyrev [KoLeSm] consider a 
three-dimensional polaron in a strong magnetic field B and with some fixed coupling a > 0. In the limit 
|B I —> oo, they argue that the dynamics of the three-dimensional magneto-polaron is equivalent to that of a 
one-dimensional strong-coupling polaron with the large coupling constant a' = (a In |B|) /2. 

But Kochetov et. al’s argument- reminiscent of Pekar’s heuristic justification [Pek] of his Produkt-Ansatz- 
has been rigorously verified only recently by R.L. Frank and L. Geisinger in [FrGs]. They calculate the exact 
ground state energy of the magneto-polaron and use the one-dimensional version of the strategy developed 
in [LiTh]. We present this one-dimensional version in full detail in order to calculate the exact ground 
state energy of our model. We also think this calculation conveniently complements Frank and Geisinger’s 
argument in Section 6 of their paper [FrGs]. 
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Theorem 1. For the ground state energy Ect{V) in 

Eo.{V) 


lim 

a—>-oo 


and the Pekar energy e(V) in P2.5\) . 
= e{V). 


The first natural task is to see that a true ground-state wave function exists. While there are existence 
results in the literature (see [Sp] [Sp-2], [Sp-3], [GrLo] and [GrLo-2]), we circumvent this important issue by 
using our result in Theorem 1 to define the notion of an approximate ground-state wave function. 


Definition 2. A wave function S .F ® L^(]R) is an approximate ground-state wave function for the 
localized Hamiltonian Ha{V) in (12.ip if H^'all.FigiL^i-R) = 1 and 

- E^{V) = o{a^). (2.7) 


If a true ground-state wave function exists, it is obviously also an approximate ground-state wave function. 
From Theorem 1 we see that Pekar’s product wave function in (12.31) is also an approximate ground-state 
wave function. We will show that any approximate ground-state wave function of our model converges in 
the weak sense to the unique minimizer of the corresponding Pekar functional given in (12.61) . 

Because the potential V{x) is symmetric decreasing and vanishes at infinity, it can easily be shown using 
the direct method in the calculus of variations that a minimizer exists for the minimization problem in (j2.5l) 
(see Theorem 8.6 and Ghapter 11 in [LiLo]; see also [Gh]). The argument is very similar to our proof of 
Lemma 5 below. A minimizer is also real-valued and positive (cf. Theorem 7.8 in [LiLo]). 

To argue uniqueness we consider the Euler-Lagrange equation. Any minimizer of the problem in (12.Sp is 
a positive, real-valued solution of the eigenvalue equation 


- u" - 2u^ -Vu = Xu, (2.8) 

where A G K is a Lagrange multiplier corresponding to the constraint u^dx = \ in (12.51) . It can be argued 
with the standard tools in regularity theory (Theorem 11.7 in [LiLo]; see also [Gh]) that the solution of (12.81) 
is in fact a classical solution. If there are two minimizers for the problem in (12.51) , then they each satisfy (12.81) 
with possibly different values of A. To show that our one-dimensional Pekar functional from (12.51) admits 
a unique minimizer, we must see that over all non-negative u G L^(IR.) with fj^u^dx = 1 there is only one 
pair (X,u) which satisfies the equation in (12.81) . We see this from Theorem 3 in [JjSt] and Theorem 2.1 in 
[MlStTr]: 

Proposition 3. IfVG C^(R) is a nonzero, nonnegative and symmetric decreasing function that vanishes at 
infinity, then the one-dimensional minimization problem in flD for the Pekar energy e{V) admits a unique 
minimizer. 

We briefly describe their proofs. First, consider the lowest eigenvalue of the linearization of the equation 
in (12.8p : 


Aq = inf 



V{x)u^dx : u G 77^ (K) and 


f 


u'^dx = 1 


(2.9) 


In 1999, H. Jeanjean and G.A. Stuart bifurcate (see Theorem 3 in [JjSt]) a unique continuous curve of 
solutions u G C'^((—oo, Aq), A) such that for each A G (—oo, Aq), the pair (A, u(A)) is a solution of (12.81) . They 
argue that all solutions of (12.81) belong to this curve: 


{(A,u(A)) : A G (—oo, Aq)} = {(p, u) : iu. i/) is a solution to ([2.81) 1. 

This means that for each —oo < A < Aq, there is a unique u{X) that satisfies (12.81) . In 2003 J.B. McLeod, 
G.A. Stuart and W.G. Troy show (Theorem 2.1 in [MlStTr]) that ||u(A)||i 2 (][{) decreases as A increases from 
—oo to Aq. Therefore, there is only one pair (A,u(A)) that satisfies (|2.8p with ||it(A )||2 = 1, and this is the 
unique minimizer of our problem in (12.5p . 

After bifurcating a curve of solutions from the lowest eigenvalue Aq in (12.91) . the strategy in [JjSt] is to 
repeatedly use the implicit function theorem at the positive solutions to get a global branch containing all 
the positive, real-valued solutions of the equation in (|2.8I) . As it was pointed out to the author by Professor 
G.A. Stuart, the uniqueness result is limited to one-dimension, because the authors in [JjSt] can only see 
how to check the hypotheses of the implicit function theorem in the one-dimensional case. Local bifurcation 
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at the lowest eigenvalue is possible in all dimensions since the eigenvalue is always simple (cf. [CrRa]), but 
global continuation is the main problem. 

With this uniqueness result we have a well-defined notion of convergence, and we now state our main 
result. 

(In this paper, we say a Borel measure W is bounded if J^W{x)dx = 1. An example is the Dirac-delta 
measure.) 

Theorem 4. Let uv be the unique minimizer of the minimization problem in and let iIIq, G T (8)L^(K.) 

be any approximate ground-state wave function of the Hamiltonian in S2.1\} . Then: 


lim 

a—>-oo 


-11^' , 
a \a 


id 


W{x)dj 


\uv\'^W {x)dx 


( 2 . 10 ) 


for any bounded Borel measure W (x) on the real line. 

2.1. Strategy. Let 5 be a real parameter and W{x) some bounded Borel measure as above. The main idea, 
suggested to the author by Professor R.L. Frank, is to perturb the Hamiltonian Ha{V) in (12.11) : 


Hc,{V-GSW)=Hc{V)-Ga^SW{ax). ( 2 . 11 ) 

Imitating the proof of Theorem 1, we calculate the exact ground state energy of the “perturbed Hamiltonian” 
in (imi) : 

EaiV+ SW) = a‘^e(y+ 6W), (2.12) 


e{V5W) =mii^v+sw{u) ■. J u^dx=l^, (2-13) 

where 

£v-\-6wiu) = Sviu) + S f W{x)\u\^dx. (2-14) 

Jr 

Of course the Pekar functional £v was already defined in (12.51) . and from Proposition 3 we know it admits 
a unique minimizer uy- 

Let be an approximate ground-state wave function of the Hamiltonian Ha(V). Since ika is not 
necessarily the ground-state wave function of the perturbed Hamiltonian Ha{V + 5W)^ a simple application 
of the variational principle yields 


E^{V + 6W) < [H^{V + 6W)] 


< Ec,{V) + o(a2) W(aa:)4'„), 

where the inner product is on E ® Lf (K). 

Dividing by and using Theorem 1, (|2.12l) and (12.131) we see that 


eiV6W) - e{V) 

S 


lim 

a—^oo 


/fV.(-) 

Jr a Va/ 


W{x) dx. 


(2.15) 


Then our main result (12.101) in Theorem 4 above follows directly from (12.151) if we can differentiate the 
“perturbed Pekar energy” e{V + SW) at <5 = 0: 


Lemma 5. Suppose the one-dimensional minimization problem in h2.5\) for the Pekar energy e(V) admits 
a unique minimizer uy- Let W(x) be any bounded Borel measure on the real line, and let e{V + 5W) be 
the perturbed Pekar energy defined in i2.13\) with some real parameter 6. Then the map 5 i-A e{V + 5W) is 
differentiable at 6 = 0 and 


jL 

Ts 


e{V + 5W) 

5=0 


/ 

JR 


W{x)\uy\'^dx. 


The proof of Lemma 5- which we provide in Section 4 below- again follows from a simple application of 
the variational principle and a compactness argument (Theorem 8.6 in [LiLo]). 
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3. Exact Ground State Energy: Proof of Theorem 1 


In this section we will calculate the exact ground state energy of our one-dimensional model using the 
strategy developed in [LiTh]. 

We use {H) to denote the expectation value of the operator on F (g)L^(]R). Likewise, |Ill'll should be read: 

The main idea in this proof is to work with a rigorously justified approximation that the electron only 
interacts with finitely many phonon modes and to then use coherent states to arrive at a lower bound that 
agrees- to the leading order in a- with Pekar’s upper bound in (|2.4|) . This already departs from the physical 
picture offered by Pekar’s ansatz, that the phonons cannot follow the electron. 

The general utility of coherent states for obtaining rigorous (lower) bounds is discussed in [LiSrYg] and 
[Li-3]. 

3.1. Ultraviolet Cutoff. We ignore large modes in the Frohlich Hamiltonian and work on a reduced mode 
space {k : \k\ < K} with a cutoff Hamiltonian 

Hk = {1- e)p^ + ^ alak - (^) " ^ - a^Viax) (3.1) 

\k\<K \k\<K 


The parameters e and K will be chosen at the very end of the computations. 

We bound the energy of the Frohlich Hamiltonian from below with the energy of the cutoff Hamiltonian 
Hk in (13.11) . We observe 

Ho^iV) =HK + ep^+ Y, ■ (z) ' ^ ( 3 . 2 ) 

|fc|>A: \k\>K 


and we will arrive at a lower bound by making an estimate on the interaction term in (13.21) . 

We use a standard commutator identity ([p, a/cC®^’^]) = k{ake'^^^) and work with two operators Z = 

{tY E ^andZt=(f)^ S 4^. 


For any 0 < e < 1, 


O’ i: ( 


ake^^^+aie-^Y )={[p,Z-Z^\) 


|fc|>Ar 


< 2(p2)= {-{Z-Z'^f)^ 

<i{t?) + i(zz' + z'z) 

= U/> + )<[z.zt]> + l(z'z) 


/ 4q; 8ck / V—^ 

^-^'Z) + 7K + 7Ki S 


a\ak 


(3.3) 


(3.4) 


\\k\>K 


Above, (13.31) is immediate from the positive definiteness of {Z -|- Z^)^. To arrive at dSH), we make the 
following estimates on (Z^Z) and ([Z, Zl]): 


kk' 
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< 


01 : 


\k\>K ' ' / 


.t. 


\|;s|>a: / 


Since [afe,a^,] = 4fc', 


([2.2']> = (i) E E - 

Ifc^iC \k'\>K 


gj(fe '=')^(afeaj^, - aj^,afc) 


fcfc' 


^(f 

We can now construct a lower bound from (j3.2ll : 

^ ^ \|fc|>Ar / ' ' 

Clearly, we require from our parameters e and ilT that 1 = |^- We now arrive at our lower bound: 

In sharp contrast to the three-dimensional computation performed in [LiTh], our error term, — i, does not 
depend on the cutoff parameter K. 


3.2. Localizing the Electron. We will bound from below the ground state energy of the cutoff Hamiltonian 

Hk- inf (ih, Hk'^) (given in ([0) 1. 
ll'flNi 

Here, (Ai?) > 0 is a parameter whose specific value will be given at the very end of the computations. 
We will denote by inf ' (di, Hk'^) the infimum taken over all wave functions whose electronic coordinate is 


localized in an interval of length 
by at most AE: 


TT 

(A£;)^ ' 


This restriction, we argue, increases the ground state energy of Hx 


inf Hx'i') > inf ' - (AE) 

111 - 11 = 1 ' ll'P||=i 

Let lilt'll = 1 and E = (4', Hx'i')- We define (j){x) = cos (^{AE)ixj on its support i 
and write (j)y{x) = 4>{x — y). To argue (13.51) . it suffices to show for some y G K, 

A direct calculation gives (((/>y'^), Hk {(fy^)) dy = J (0)^ -I- Ecjy^dx and 


[ {{(cfy^), HK{f>y'it))-{E + AE){cfy^, cfy^))dy 

Jr 


(3.5) 

TT TT \ 

2{AE)i ’ 2(AE)5 J 

(3.6) 


= (0')2 _ {AE)cj)'^dx = 0 

J - 

2(AS) 2 

since (AE) is the Dirichlet energy of (f>. So there exists some y G M such that (|4.6ll holds. From now on, we 
consider the electron to be localized in some interval of length —. 

(AE)2 
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3.3. Block Hamiltonian. We now decompose our finite mode space into finitely many blocks: {k : |A:| < 
if} = each block Bn contains “PL” modes where “ max \ki — %| = P” is the size of each block. 

n ki^kj^Bn 

On each block P„ we analogously define block annihlation and creation operators: Ab^ = —X) ^k 
and ylg = —3 -t- ^ a\. Clearly, [Ab„^,^b ] = ^rnn- On each block P„, we see from the Cauchy-Schwarz 
inequality that 

^ aj^Ofe. (3-7) 

feeSn 

For reasons that will become clear in the next stage of the computation, we now want to work with the 
approximation that the electron interacts with atmost one mode in each block P„. For this approxi¬ 
mation to work, we make the following estimate on the interaction term of our cutoff Hamiltonian: for any 
parameter 0 < 5 < 1 and any mode in each block Bn, completing the square yields 




Bn k€Bn 


\ \k\<K 




iksn ^ ikx I ^ 


Bn kGBn 


< ( <^ II + (f) ^ H H 1^ - 


P|a:P 


\k\<K 


\ \k\<K 


B„ keB„ 

2aKP^'K^ 


S{AE) 


(3.8) 


To arrive at (13.8p . we used the rigorously justified approximation (see p.Sp ) that the electronic co-ordinate 
is localized in an interval of length —, 

(AE)2 

The parameters P and 0 < 5 < 1 will be chosen at the very end of the computations; the specific mode 
kB„ in each block will be chosen in the next stage of the computation. 

Now we bound the ground state energy of the cutoff Hamiltonian (with the condition that the electronic 
co-ordinate of the ground-state wave function is localized) inf ' (dt, idifT), from below, using the energy 

of the block Hamiltonian: 

= (1 - e)p2 -k (1 - <5) ^ A^Ab^ - {Pa)i “ o:^V(ax). (3.9) 


Bn 


Clearly, 


Hk = ( (1 - e)L^ + H alak 

\k\<K 


(i)'i:i:( 




Bn k&Bn 


^V{ax) + (I) ' H H [ak ^ 


Bn k^Bn 


^ o-Block^r,- 2aifp27r2 

>H^ {{kBn))-^^^, 


(3.10) 


in the sense of expectation values. To arrive at (13.101) we simply used the estimates from (13.71) and (13.8F 
We summarize: 


inf' y, Hk^)> inf' sup /d-, Llf ({fcsj) d' 


2aKP^Tr^ 


S{AE) 












GROUND-STATE WAVE FUNCTION OF THE ONE-DIMENSIONAL POLARON 


10 


3.4. Coherent States. We now work with the block Hamiltonian ({^b„}) from (13.91) and the block 

creation and annihlation operators constructed in the previous stage of the computation. For each block Bn 
we define a block coherent state indexed by some zb^ € C: 

NsJ = |0 bJ , 

where |0 _b„) denotes the vacuum state in block Bn, i.e., Ab^, |0b„) = 0. We write 

1^) 

Bn 

a tensor product of the coherent states corresponding to each block Bn- 

One can verify that for each block Bn, the coherent state |z_b„) is the eigenstate of the corresponding 
block annihlation operator: 

\ZBn) = ZBn \ZBn) ' 

The commutator relation [Ab„^, Ag^] = 6mn, together with the resolution of identity formula 

1= k) {z\Y[dzBn,(t^, 

Bm 

yield the convenient representations: 

ABn = / ZBn\z) {z\Y\_dZBn,dzE^ 

Bm 


Ai A 


Bn^Bn 


\ZBn 


- l) k) {z\Y[dzB„,d:^ 


Denoting 4'z(x) = {z \ 'I')pjionon(f^® inner product only in the phonon coordinates), we recast the energy 
of the block Hamiltonian ({^s„}) in the following form: 


4-, 


{{kBn})^)= [ ('k., WeH dzBn^d^ 

D 


(3.11) 


where (•, •)e is the inner product only over the electronic coordinate, and is the Schrodinger operator: 

K = {1- e)p^ + ^ [(1 - <5) {\zBn I' - l) - (Pa)i - a^V{ax) (3.12) 


Since 


4'z)b(1 - S)z^ - (Pa)5 (4-,, {^b^ 


{Pa) 2 


(1 -5)(4'^,4'^) 


('kz, 


— ikSn ^ 


'^z)e >0, 


completing the square yields 

(l-5)(^.,^,)£;kB„p-^(Pa)5(^z,e-*'=®""4',)B-ZB„(Pa)^(4'z,P'=®""^.)£; 


> 


-(Pa)|(4'.,e-*'=^""T,)B| 


(1-(5)(4'^,4'^)b 

The advantage of constructing a Block Hamiltonian in the previous subsection is that the energy error we 
incur for disregarding the “-1” term in (13.1211 is proportional to the number of blocks: a finite value. 

In the following calculation, in each block Bn we choose a mode /Cb„ such that 


\{^z,e = min |(4'z,e *'="^4',,) 


k^Bn 


E\ 


As seen in (13.151) below, we also make use of the continuum approximation L J dk permitted by our 

model. We now proceed to extract the Pekar energy functional from (j3.11|l : 


sup ({fcsU) ^ 
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> 


('I'*. 


_ (1 _ 

/(*., »dE X (^(1 - 


^)eI 


El 


(3.13) 


2K 

P 


a 1 

(1-5) L ^ 






El 


, {'F,, V(ax)^,)^ 

('I'.. 'I'.>g 


Ub^ dzB^d^- (l-<5) — 


> 


/(». *.>Ex((l-e,<;;^-^iE^ 


P 

— ikx 


'I'.) 


'I'-)eI 

- 2 - 


2 y{oix)^^ 

““ ('F.. 'I'.)g 


El 


n_B„ dzB^dzEZ- (1 - <5)^ 







\ 

('F^, 'F^)g 

n 


dfc-l- 


2 ('F., V{ax)^.)y^ 
(■F^, 'F^)g 


ris^ dzBmd ^ z ^- { 1 - 6 )^ 

'F^l^ 




> 


= /(U, «dEX 

Hr™ dzB^dz^- (1 - 6) 

f (^z, ^z)e X (1 - e) (-T^TTny— “ (i_5“(l-e) / 


(i2;+ 


(3.14) 


(3.15) 


(3.16) 


E 

2K 

P 


<3'*. 'F,) 


dx+ 


El 


y(aa;)'F^) 


E 


> 


(l-e)"(l-5)" ('I'.. 'I'.>E 

a^e(E) 


Ob dzB^ctz^ - {I - 5) 


2K 


J i'i’z, «'z)e 


(l-e)(l-5)2 

a'^e(y) 


2K 


Y[dzB^dzBm - (1 - d) — 


B„ 


(3.17) 


(l-e)(l-<5)2 ^ ’P^ 


where e(E) is the Pekar energy, defined in (12.51) . Above (13.161) follows from the Plancherel’s theorem, and 
(13.171) follows from the scaling properties of our Pekar functional. 

3.5. Controlling the Error Terms. We have the following upper and lower bounds on the ground-state 
energy: 


a^eiV) > E^{V) > 


e{V)a^ 2K 1 2aKP^n^ 

(l-e)(l-5)2 - (1 - - 2 - - S{AE) 


= a^e{V) — 


a 


2x2 n„.2x SA I 16a^5 8a^5^ ' 

0 za 0 — -I — —^ 


1-26 + 6'^-^ 


8a _i_ 1 6a(5 8a(5^ 

K K 


error-term 


since we noted, when using an ultraviolet cutoff above, that the parameters e and K must satisfy the coupling 
relation: e = ^. We now choose specific values (in orders of a) for the parameters 77, (5, P and AE so that 
the error-term above is of an order less than c?^ while satisfying the following constraints: 0 < 5 < 1 and 
P < K when a 1. In an attempt to make the error-term as small as possible, we have chosen 


6 = cia ,K = C 2 a« , P = caa^ and AE = C4Q;<‘9 . 


From (13.181) we conclude 
This proves Theorem 1. 


Q^e{V) > Ea{V) > a^e(V) - Ca^. 
















































GROUND-STATE WAVE FUNCTION OF THE ONE-DIMENSIONAL POLARON 


12 


4. Differentiating the Pekar Energy: Proof of Lemma 5 
Proof. Given any e > 0 and for any finite parameter S > 0, choose us € so that ||m 5||2 = 1, and 

[ {\us\^-v{x)\usf-\us\^ + SW{x)\usf)dx-e{V + SW) <e (4.1) 

Jr 

Since uy is not necessarily the minimizer corresponding to the perturbed Pekar energy e{V + 5W)^ 

[ {\u's\^ - V{x)\us\'^ - \us\^ + SW{x)\us\‘^) dx < f (|w'yP “ Vix)\uvf - + SW{x)\uv\‘^) dx, (4.2) 

Jr Jr 


and 


^ I \uv\^W{x)dx + e{V) > e{y + 5W). 

Jr 

Since us is not necessarily the minimizer corresponding to the Pekar energy e{V), 


(4.3) 


f {\Vus\^+ V{x)\us\^ — \usf + 6W{x)\us\^) dx > e{V) + S f \us\‘^W{x) dx, (4.4) 

Jr Jr 


and 


j luvl"^W{x)dx > 

Jr 


e{V + 6W) - e{V) 


> 


j \u5\'^W{x)dz 

Jr 


(4.5) 


In parallel had we chosen (5 < 0 we would have arrived at the relation in (|4.5I) with the inequalities only 
reversed. 

It suffices to show 

(4.6) 


lim f \us\'^W(x) dx = f \uv\'^W{x) dx, 
^-*•0 Jr Jr 


where W (x) is a bounded measure on the real line. We show this by arguing that us converges uniformly to 
uv on the real line. 

In the limit 5^0, the set of functions {rt^} from (j4.1|) above is in fact a minimizing sequence of the Pekar 
functional Ev (see (j2.6ll l with the Pekar energy e{V). For 5 small enough, since e{V) < 0, 


1 > 1 


i{V)>£v{us) = J (^\us\'^ - lusl'^ - V{x)\us\'^'^ dx 


>\\\ns\\l-\- 


(4.7) 

(4.8) 


In (j4.7l) we use that P € (R) and is a symmetric decreasing function. In (j4.8ll we use the one-dimensional 

Sobolev inequality, whence we have the uniform bound 

Il'*^(5ll2 < 3 + 2||P||oo. 

The sequence {rti} is thus uniformly bounded on and {rti}- or a subsequence thereof- has a weak 

limit u £ 

Since P is a symmetric decreasing function, using rearrangement inequalites we see that given any e > 0, 
there is a compact set Tfg and a parameter > 0 such that 

{a; £ R : > e for all 0 < i5 < Dg} C (4.9) 

Since P vanishes at infinity, we see that given any e > 0 there is a compact set Tfe such that 

{a: £ R : P(a:) > e} C K^. 

The compact set obviously has finite measure, so appealing to Theorem 8.6 in [LiLo] we conclude that the 
sequence {rt^} converges to u strongly in L^^Kf) for all p < oo, and that {u^} converges pointwise uniformly 
to u on K^. The latter convergence result follows from the fundamental theorem of calculus. 

Thus, 

ix — f Irtl^da;. 


■ f \us\'^V{x)dx — ( \us\^dx —)• — f \u\'^V{x)da 

Jr Jr Jr 


This establishes the weak lower semicontinuity of the Pekar functional Ey, and the weak limit u is in fact the 
minimizer of Ev- But we see from Proposition 3 above that Ey admits a unique minimizer uy, so u = uy. 
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Thus Us converges pointwise uniformly to uy on and from (14.91) we easily see that the convergence can 
be made uniform on the entire real line. 

The relation in (14.6p now follows for any bounded measure W(x) on the real line, and the perturbed Pekar 
energy e(V + SW) is indeed differentiable at 5 = 0: 


e(V + 6W) - e(V) 
hm-;- 

il-i-O 0 


\uv\‘^W{x)dx. 


□ 


5. Proof of Theorem 4 

In the proof, (•, •) should be read as the inner product on 0 L^(R). 

Proof. Recall the perturbed Hamiltonian Ha{V + SW) defined in (12.111) : 

HaiV + SW) = HaiV) + a^SWiax) 

and its ground state energy 

Ea{V + SW) = a^e(y + SW) 

discussed in (12.121) - (12.141) . 

Let Tq be an approximate ground-state wave function (see Definition 2 above) of the Hamiltonian HaiV) 
given in dUD. From the variational principle. 


so that 

and 


E^iV + SW)-{^c.. H^iV)^c.) <a^S{^c., Wiax)'^^), 
a-^E^iV + SW)-{^^, VF(ax)T„), 


liminf a {Ea{V + SW) — EaiV)} < S bP(aa;)'I'Q,)') . 

Oi—^OO \ Cl—¥00 J 

From Theorem 1 and ()2.12p above, 

e(P -I- (HP) — e(H) < (5 ( liminf ( / W{x) — H^'a f— ) 

V \Jr C( \a/ 


dx 


For S > 0 


and 


e{V + SW) - e{V) 


lim 

(5—>-0+ 


e{V + SW) - e{V) 


< liminf f Wix) ( —llTo, f—) ll^r ) dx, 

Jr V« / 

< liminf ( [ W (x) — H^'a f—) ||^ dx). 
a^-oo \J^ \_a \a/ \ J 


From Lemma 5, 

f W(x) Iwypda; < liminf ( [ W(x) 

Jr a->-oo 

For d < 0, we similarly arrive at the following relation: 

I lT(a:)|Mypda; > limsup ( f W{x) 


a va> 


o 

G) 


We thus conclude: 


[ W{x)\uvfdx = lim ( f W{x) 
Jr q^oo V Jg 


a \aj 


a a 


dx I . 


dx 


dx 


which is our desired convergence relation for any approximate ground-state wave function with W a bounded 
measure on the real line. □ 
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